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Abstract 

We present an old and regretfully forgotten method by Jacobi which allows one to 
find many Lagrangians of simple classical models and also of nonconservative systems. 
' r*j | We underline that the knowledge of Lie symmetries generates Jacobi last multipliers and 

Q_f each of the latter yields a Lagrangian. Then it is shown that Noether's theorem can 

identify among those Lagrangians the physical Lagrangian(s) that will successfully lead 
to quantization. The preservation of the Noether symmetries as Lie symmetries of the 
corresponding Schrodinger equation is the key that takes classical mechanics into quantum 
mechanics. Some examples are presented. 
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1 Introduction 

The inverse problem of calculus of variation has attracted a lot of interest since in the second 
half of the XVIII century Euler [8] and then Lagrange [IT] introduced the direct problem, 
namely the idea of linking the solution of a differential equation to the maximum/minimum 
of a functional, the celebrated problem of the brachistochrone being indeed the most famous 
classical example. It will take hundred of pages to cite all the papers and books that have 
been published since up to date. Most authors mark the birthdate of the inverse problem with 
the 1887-papers by either Helmholtz [11] or Volterra [33]. Some other especially among the 
Russian speaking researchers pushes the date slightly back to the 1886-paper by Sonin [32J. 
Very few recognize the seminal work by Jacobi, namely his 1845-paper [14] and his 1842-1843 
Dynamics Lectures published posthumously in 1884 [15], available in English since 2009 [16J, 
where he links his last multiplier to the Lagrangian for any even-order ordinary differential 
equation (ODE). Actually both Volterra and Sonin recognize the contribution of Jacobi last 
multiplier in their papers, Sonin more explicitly than Volterra since he showed that his own 
method involves the Jacobi last multiplier (p. 10 in [32]). 
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The method of Jacobi last multiplier was enhanced when Lie determined the link with his 
symmetries [18J, a link very easy to implement that allows to derive many multipliers and 
therefore Lagrangians. 

It is known that a single second-order ODE admits different many Lagrangians |34| . but 
so far there is not a method that can discern the physical Lagrangian among them, although 
some were proposed, e.g. |6j. The same is true for systems of second-order ODEs that admit 
more that one Lagrangian. It has been shown that some systems of second-order ODEs do 
not admit a Lagrangian [7J although in [27] the Lagrangian of some of those systems were 
determined by following Bateman's statement [2j, namely finding a set of equations equal in 
number to a given set, compatible with it and derivable from a variational principle without 
recourse to any additional set of equations. 

We propose that the physical Lagrangian should be the one that admits the highest possible 
number of Noether symmetries |26j , [23] , |24j . It was proven in |10] that the maximal dimension 
of the Lie symmetry algebra of a system of n ODEs of second order is n 2 + 4n + 3, and the 
highest number of corresponding Noether symmetries is (n 2 + 3n + 6)/2. In particular for a 
single second-order ODE the highest number of Noether symmetries is five |20| . 

Consequently we conjecture that the passage from a classical system to its quantum ana- 
logue should preserve exactly those Noether symmetries, namely the Noether symmetries of 
the physical Lagrangian shall become the Lie symmetries of the corresponding Schrodinger 
equation [26]. [23]. [24]. 

In this paper after recalling the method of Jacobi last multiplier and its link to Lie sym- 
metries, we present the simple example of the one-dimensional free particle: ten inequivalen10 
Lagrangians are presented, their Noether symmetries identify two Lagrangians that admit the 
highest number of Noether symmetry: one independent and one dependent on time. Then the 
corresponding Schrodinger equations are obtained. Also the case of the nonlinear second-order 
Riccati equation is illustrated and its corresponding Schrodinger equation obtained. 



2 Jacobi last multiplier 

The method of the Jacobi Last Multiplier |12j . [13] . [14] . [15] provides a means to determine 
all the solutions of the partial differential equation 

n df 
Af = ^2a i (x 1 ,...,x n )— =0 (1) 

i=l ' 

or its equivalent associated Lagrange system 

ffl _ ^2 _ _ ffn ^) 

In fact, if one knows the Jacobi Last Multiplier and all but one of the solutions, then the last 
solution can be obtained by a quadrature. The Jacobi Last Multiplier M is given by 

d(xi,x 2 , ...,x n ) 



1 Namely they do not differ by a total derivative. 
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where 



d(f,UJ 1 ,UJ 2 ,. . . ,UJ n ^i) 

d(x 1 ,x 2 ,...,x n ) 



det 



df 


of 


dx\ 
dux 


dx n 
dbJi 


dx\ 


dx n 


du) n -\ 


du n -i 


. dx\ 


dx n . 



(4) 



and ui, . . . ,u) n —i are n — 1 solutions of (HJ or, equivalently, first integrals of ([2]) independent 
of each other. This means that M is a function of the variables (x±, . . . ,x n ) and depends on 
the chosen n — 1 solutions, in the sense that it varies as they vary. The essential properties of 
the Jacobi Last Multiplier are: 

(a) If one selects a different set of n — 1 independent solutions r/i,... ,r/ n -i of equation ([T|), 
then the corresponding Last Multiplier N is linked to M by the relationship: 

N = M _d( m ,...,T h - 1 ) 



d(uji,...,u) n -i)' 
(b) Given a non-singular transformation of variables 

t : (xi, X2, • • • , x n ) > (x^, x 2 , • • • , ^ n )j 
then the Last Multiplier M' of A'F = is given by: 

r d(xi,x 2 , ■■■ ,x n ) 



M' = M 



where M obviously comes from the n — 1 solutions of AF = which correspond to those 
chosen for A'F = through the inverse transformation t _1 . 

(c) One can prove that each multiplier M is a solution of the following linear partial differ- 
ential equation: 

d(M ai 



E 

8=1 



dxi 



0. 



or equivalently: 



(5) 



(6) 



viceversa every solution M of this equation is a Jacobi Last Multiplier. 



(d) If one knows two Jacobi Last Multipliers Mi and M 2 of equation ([TJ) , then their ratio is a 
solution uj of dl]), or, equivalently, a first integral of ([2]). Naturally the ratio may be quite 
trivial, namely a constant. Viceversa the product of a multiplier M\ times any solution 
uj yields another last multiplier M 2 = M\uj. 
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Since the existence of a solution/first integral is consequent upon the existence of symmetry, an 
alternate formulation in terms of symmetries was provided by Lie [18] , [19] . A clear treatment 
of the formulation in terms of solutions/first integrals and symmetries is given by Bianchi [3]. 
If we know n—1 symmetries of ([1]) / ([2]) , say 



r i — ^ (%x j • • • 3 x n ) X j , i — 1 , ti 1 , 

3=1 



(7) 



Jacobi's last multiplier is given by M = A , provided that A ^ 0, where 



A = det 



ii,i 

£n-l,l 



£l,n 



in- 



l,n 



(8) 



There is an obvious corollary to the results of Jacobi mentioned above. In the case that there 
exists a constant multiplier, the determinant is a first integral. This result is potentially very 
useful in the search for first integrals of systems of ordinary differential equations. In particular, 
if each component of the vector field of the equation of motion is missing the variable associated 
with that component, i.e., dai/dxi = 0, the last multiplier is a constant, and any other Jacobi 
Last Multiplier is a first integral. 

Another property of the Jacobi Last Multiplier is its (almost forgotten) relationship with 
the Lagrangian, L = L(t,q,q), for any second-order equation 



i.e. [15] (Lecture 100, [34J 



M 



d 2 L 
dq 2 

where M = M(t, q, q) satisfies the following equation 



0. 



Then equation ([9]) becomes the Euler-Lagrange equation: 



d fdL 



dt \dq ) dq 



8L 



(9) 
(10) 

(11) 

(12) 



The proof is based on taking the derivative of (|12p with respect to q and showing that this 
yields (jlip . If one knows a Jacobi last multiplier, then L can be easily obtained by a double 
integration, i.e.: 

L = J (J Md^j dq+f 1 (t,q)q + f 2 (t,q), (13) 



2 Jacobi's Lectures on Dynamics are finally available in English [16] . 
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where fa and fa are functions of t and q which have to satisfy a single partial differential 
equation related to ([9]) [25]. As it was shown in [25], fa, fa are related to the gauge function 
g = g(t,q). In fact, we may assume 



h 
fa 



Oq 

| + A(M) 



(14) 



where fa has to satisfy the mentioned partial differential equation and g is obviously arbitrary. 
We remark the importance of the gauge function in order to apply Noether's theorem [21] 
correctly. Therefore we do not annihilate the gauge function. 

In [25] it was shown that if one knows several (at least two) Lie symmetries of the second- 
order differential equation ([9]), i.e. 



Tj = Vj(t, q)8t + Gj(t, q)d q , j = 1, r, 
then many Jacobi Last Multipliers could be derived by means of ( 

1 q 



i.e. 



1 



det 



Vn G r , 
Vm. G 



F(t,q,q) 
G V 

t q t 



g; 



t 1 t 

with (n,m = 1,7"), and therefore many Lagrangians can be obtained by means of (|13|) . 



(15) 



(16) 



3 Lagrangians for the free particle 

It is well-known since Lie's seminal work [19] that the equation of a free particle 

q = (17) 

admits an eight-dimensional Lie symmetry algebra^, sl(3,M), generated by the following op- 
erators: 



Xi = qtd t + q 2 d q , X 2 = qd t , X 3 = t 2 d t + qtd q , X 4 = qd q , 



X, = tdt 



x 6 = d t , 



X 7 = td, 



9! 



^8 = Og. 



(18) 



An obvious Jacobi Last multiplier (JLM) of (11 7f) is a constant since q does not appear in its 
right-hand side. This also implies that any JLM is a first integral of (I17D . We can found this 
trivial JLM and many others by calculating the inverse of the determinant of the matrix (|16|) 
for all the possible combinations of two different operators in (|18p . It results that ten different 



3 Also reported in [5]. 
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JLM and consequently as many Lagrangians, by means of (jl3[) . can be obtainecQ, i.e. 



M 15 = -— — => L 15 = ±(log(tq-q)-log(q)) + -^(t,q) 



q 2 ™"* " ' dt 

Mir, - ~ - L l{ , ( 4 - (log(g) - log(ig - q)) + -£(t, . 

M i8 = -— - L 18 = -llog{q)-(~-l)log(tq-q) + ^(l+log(q)) + -£(t,q) 

qitq-q) q \t qj t dt 

1 r 1 d^/ N 



(*<?- 


g) 3 


1 




g(ig - 




1 




em- 


-9) 








g) 2 


1 






-9) 



4 => L 28 = -log(g) + ^( 
g^ dt 



M -=i^ - ^=(i-i)(iog W - 9 )-i )+ ^( t , 9 ) 

M 48 = -i =► L 48 = g(l-log(g)) + ^j(t,g) 

M 87 = l L 87 = i g 2 + ^|(t,g) (19) 

The ten Lagrangians are NOT linked by the gauge function g = g(t, q). We note that seven of 
the matrices (|16p have determinant equal to zero, i.e.: 

Ai 2 , A14, A 24 , A 35 , A 37 , A 57 , A 68 . (20) 

This is in agreement with the application of Jacobi's method to the linear harmonic oscillator 
|25j that yielded fourteen different Lagrangians and three matrices having determinant equal 
to zero. 

Applying Noether's theorem yields that both L13 and L 87 admit five Noether point symme- 
tries, the maximum possible |20j in the case of equation (117p . The main difference between 
the two Lagrangians is that one Lagrangian only is independent on time and is the traditional 
Lagrangian, namely the kinetic energy L 87 . Moreover the two Lagrangians admit different 
Noether point symmetries. In fact Lagrangian L13 admits the following five Noether symme- 
tries and corresponding first integrals of equation (|17j) 

Xi => InU 



q-tq 



4 Since M nm — ~M mn we have arbitrarily chosen the sign as we wished. 
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while Lagrangian L^-j admits the following five Noether symmetries and corresponding first 
integrals of equation (fT7|) 



x 3 = 


=> In 3 = 


-3(9- 


tq? 


X4 + 2X5 = 


=> In& = 


-q{q - 


tq), 


x 6 = 


=> In 6 = 


2 q > 




x 7 = 


=4> Inj = 


g - tq, 




X s = 


=> In 8 = 


-q- 





(22) 

4 Schrodinger equations for the free particle 

The Schrodinger equation for the free particle is 

2#t + 1>xx = 0. (23) 
We show that this equation can be obtained by considering a generic linear parabolic equation 

2i^ t + fx(x)ip xx + f 2 (x)il> x + h{x)^ = (24) 

with fk,(k = 1,3) functions of x to be determined in such a way that equation (f2l"|) admits 
the following five Lie symmetries^] 

X 3 => = t 2 d t + xtd x + wic^, 

X 4 + 2X5 n 2 = 2td t + xd x + oj 2 d^, 

X 6 => Q 3 = d t +to 3 d^, 

Xj Q4 = td x + uj^d^,, 

X 8 => f2 5 = ^ + W5^. (25) 

where Wj = Ui(t,x,if)), (i = 1,5) are functions of t,x,ip that have to be determined. Equation 
(|24f) also admits the following two symmetries 

^6 = = a(t, a:)<9</, (26) 



3 We have identified g with x. 



with a any solution of equation (|24p itself, since any linear partial differential equation possesses 
these two symmetries. 

Using the interactive REDUCE programs [22], we obtain that 

/i = l, / 2 = / 3 = 0, (27) 

and 

= \{ix 2 — t)lf}, U>2 = ^3 = ^5 = 0, u>4 = ixip. (28) 

Therefore the Noether symmetries admitted by the Lagrangian L§7, namely the physical La- 
grangian for the free particle, yield the right Schrodinger's equation (j23|) and thus the correct 
quantization procedure is achieved. 

We now show that a quantum-correct Schrodinger equation can be derived even from the 
five Noether symmetries admitted by the time-dependent Lagrangian L\ 3 . We consider a 
generic linear partial differential equation 

fii(t,x)ip u + fi2(t,x)ip tx + / 2 2(* 5 x)ip xx + fi(t,x)i> t + f2{t,x)tp x + f Q (t,x)ip = (29) 



x 1 = 




x 2 = 




x 3 = 


» w 3 


x A - x 5 = 




x 7 = 


> w 5 



with f rs , f r ,(r,s = 1, 2), /o functions of t,x to be determined in such a way that equation ([29]) 
admits the following five Lie symmetrieqj 

xtd t + x 2 d x + widip, 
xd t + w 2 d, l i ) , 
t 2 d t + xtd x + w 3 d^, 
-td t + xd x + 104^, 

td x + w 5 d^,. (30) 

where Wi = Wi(t,x,ijj), (i = 1,5) are functions of t,x,ip that have to be determined. Equation 
(I29p also admits the two symmetries 

W 6 = ^, Wp = /3(t,x)d^ (31) 

with (5 any solution of equation (|29l) itself. 

Using the interactive REDUCE programs [22], we obtain that equation (j29|) becomes 

At 2 ^ tt + Stx^ tx + Ax 2 i> xx + 12tipt + I2x^ x + 3V = (32) 

with 

Wi = -±a^, PU 2 = W 4 = W 5 = 0, PU 3 = -\tijj. (33) 

Equation (132p is parabolic and therefore can be put into its normal form since its characteristic 
coordinates are 

£ = |, x = x, ^ = 4>{^x). (34) 



°We have identified q with x. 



s 



Thus equation (j32[) transforms into 

Ax 2 (j) xx + I2x4> x + 30 = (35) 

with solution: 

<P(t x) = ariOx- 1 / 2 + a 2 (Ox- 3 / 2 , (36) 

where a%, a 2 are arbitrary functions of £. This is obviously in agreement with the quantization 
of the free particle. 

We remark that there is some freedom when imposing equation (I29p to admit the Lie 
symmetries (|30p but in any possible cases a parabolic equation with characteristic coordinates 
(|34p is always obtained. For example another possibility is 



At 4 ip tt + 8t 2 xi/j t x + 4t 2 x 2 ip xx + At 2 (3t + x)ip t + 4tx(3t + x)ip x + (3t 2 + 4te + x 2 )^ = (37) 
with 



W 



| (l + f ) ^2 = ^l0g(t)V, W 3 = -i(t + x)0, 



W 4 = - log(t)V, W 5 = -i log(i)^. (38) 

The normal form of equation (|37p is 

4x 2 (t> xx + (12 + 4£)x</> x + (3 + 4£ + £ 2 )0 = (39) 

with solution: 

^i) = MOx- 3 / 2 ^ 2 + a 2 (Ox- 1/2 -* /2 , (40) 
also in agreement with the quantization of the free particle. 

5 Schrodinger equation for the second-order Riccati equation 

In [25] it was shown that the linearizable second-order Riccati equation, a member of the 
Riccati-chain |T], i.e. 

x + 3xx + x 3 = (41) 

possesses many JLM and therefore Lagrangians. In particular the following time-independent 
Lagrangiarjll 

Lagr = - — ± (42) 
2(x + x z ) 

was shown to admit five Noether point symmetries, i.e. 

r 2 -r 8 , r 3 -^r 7 , r 4 , r 5 , r 6 (43) 



7 This Lagrangian has also been studied in [5]. 
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among the following eight Lie symmetries admitted by equation (|41H : 



ri 


= t 3 {tx 


- 2)d t - t(xt - 


2){x 2 t 2 + 2-2xt)d x 


r 2 


= xt 3 d t 


- (xt-l)(x 2 t 2 


+ 4 - 2xt)d x 


r 3 


= xt 2 d t 


- x(x 2 t 2 + 2- 


2xt)d x 


r 4 


= xtdt - 


- x 2 (xt - l)d x 




r 5 


= xd t - 


x 3 d x 




r 6 


= d t 






r- 


= td t - 


xd x 




r 8 


= t 2 d t - 


- 2{xt - 1)8 X . 





(44) 



Lagrangian (142j) was obtained from the inverse of the determinant of the matrix ()16f) with the 
two symmetries T$ and Tq and the application of formula (|13p . 

It interesting to show that equation (|41h is linked to the following cubic-power dissipative 
equation 

d 2 x fdx\ 3 , d . 
+ 37 =° ( 45 ) 



dt 2 \dtj 

by means of the canonical form of the two-dimensional Lie algebra generated by Ts and Tq . In 
fact this Lie algebra corresponds to Type I in Lie's classification of the real two-dimensional Lie 
algebras in the plane [19], [3], namely it is abelian and transitiv^f] and therefore its canonical 
coordinates are 

d h a s . (46) 

Therefore the following identification 

r 5 = d h r 6 = d x (47) 

yields the transformation 

1 tx — 1 , , . 

t = ^ X = — (48) 

that takes (JUJ) into ([35]) . 

Equation (|45h is also linearizable and admits an eight-dimensional Lie symmetry algebra. 

Work is in progress to quantize n-power dissipative equations with the method illustrated here 
and to compare it to other proposed method, e.g. [29], [30], [3T] . 



We now quantize the second-order Riccati equation (|4Tj) by imposing equation (129]) to 
admit the following five Lie symmetries 



r 2 -r 8 = 


> Ai 


= {xt - 


- l)t 2 d t - (xt-l){x 2 t 2 


- 2 - 2xt)d x + Ai^ 


r 3 - ~r 7 = 

6 3 ' 


> A 2 


= (xt 


-l)td t -x(x 2 t 2 + l 


- 2x t) d x + \ 2 d^, 


r 4 = 


> A 3 


= xtd t 


- x 2 (xt - l)d x + A 3 9^, 




r 5 = 


» A 4 


= xd t - 


- x 3 d x + A 4 5^, 




r 6 = 


» A 5 


= d t + 







(49) 



In [28] a missprint led to the erroneously statement that and Fe generate an intransitive Lie algebra. 
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where Aj = \i(t, x, = 1, 5) are functions of t, x, ip that have to be determined. We remind 
that equation (|29p also admits the two symmetries (|31|) . 

Using the interactive REDUCE programs [22], we obtain that equation (I29p becomes 

4^t ~ 8x 2 ^ tx + 4xVx* + 8x 3 ^ - 3x 2 ^ = (50) 

with 

(tx - l) 3 



2x 



-if), A2 = —\(tx — l) 2 ip, A3 = — \{tx — l)xif) 

fx — 1 

A 4 = -ixV, A 5 = V- (51) 



Equation (|5U|) is parabolic and therefore can be put into its normal form since its characteristic 
coordinates are 

£ = i , x = x, ip = (f>(g,x). (52) 

x 

Thus equation (|50p transforms into 

^x 2 <p xx + 8x(fe - 30 = (53) 

with solution: 

x) = /3i(g)x^ 2 + /3 2 (^~ 3/2 , (54) 



where /3i,/02 are arbitrary functions of Is this the correct quantization of equation (14ip ? It 
seems so but further insight is needed especially from the experimentalists. 

6 Final remarks 



In [23] a method was proposed to overcome the deadlock of nonlinear canonical transformations 
when quantizing with the known procedures [26] . It consists of the following steps to be applied 
to a classical Lagrangian equation^: 

• Find the Lie symmetries of the Lagrangian equation 

T = W (t,x)d t + W 1 (t,x)d x 

• Among the Lie symmetries find the Noether symmetries admitted by the given La- 
grangian 

T = v {t,x)d t + v 1 (t,x)d x , rcT 

• Construct the Schrodinger equation admitting these Noether symmetries as Lie symme- 
tries 

2iipt + fi(x)ip xx + f2{x)ip x + h{x)i> = 
n = V (t, x)d t + Vi (t, x)d x + G(t, x, z/Ofy, 

• Summarizing: quantize preserving the Noether symmetries 



The physical Lagrangian turns out to admit the highest number of Noether point symmetries. 
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In [23] this method has been applied to classical Lagrangian systems. In particular it led to the 
Schrodinger equation of a known completely integrable and solvable many-body problem, the 
so-called 'goldfish' [3]. In the case of Lagrangian systems the method consists of the following 
step© 

• Find the Lie symmetries of the Lagrangian system 

N 

T = W(t,x)d t + Y / W k (t,x)d Xk 

k=l 

• Among the Lie symmetries find the Noether symmetries admitted by the given La- 
grangian 

N 

r = v(t,x)d t + ^v k (t,x)d Xk , rcT 

k=l 

• Construct the Schrodinger equation admitting these Noether symmetries as Lie symme- 
tries 

N N 

2i^ t + ^2 fkj(x)i/) XjXh + h k{x)ip Xk + h{x)^ = 

k,j=l k=l 
N 

n = v(t,x)d t + £ v k (t,x)d Xk + G(t,x, vo<^ 

k=l 

• Summarizing: quantize preserving the Noether symmetries 

In this paper we have proposed a method that takes classical equations into the quantum 
realm by means of the Schrodinger equation regardless of the number of Lagrangians that may 
classically exist. The method consists of the following steps: 

• Find the Lie symmetries of the given equation 

T = W (t,x)d t + W 1 (t,x)d x 

• Among many Lagrangians identify the ona[^] that admits the highest number of Noether 
symmetries among the Lie symmetries 

r = v (t,x)d t + v 1 (t,x)d x , rcT 

• Construct a linear partial differential equation admitting these Noether symmetries as 
Lie symmetries 

fn(t,x)ip tt + fi2(t,x)iptx + f22(t,x)ip xx + fi(t,x)rpt + f2(t,x)ip x + f (t,x)ip = 

n = V (t, x)d t + Vi (t, x)d x + G(t, x, z/))^ 

• Since this equation is parabolic, determine its canonical coordinates in order to transform 
it into its normal form, namely the Schrodinger equation 

• Summarizing: quantize preserving the Noether symmetries. 

10 The physical Lagrangian turns out to admit the highest number of Noether point symmetries. 
n They could be more than one as we have shown in the present paper. 
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